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Abstract 

In this article we present a local hidden variables model for all experi- 
ments involving photon pairs produced in parametric down conversion, based 
on the Wigner representation of the radiation field. A modification of the 
standard quantum theory of detection is made in order to give a local realist 
explanation of the counting rates in photodetectors. This model involves the 
existence of a real zeropoint field, such that the vacumm level of radiation 
lies below the threshold of the detectors. 
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1 Introduction 



Following Bell, a local hidden variables model (LHV) exists for an Einstein- 
Podilsky- Rosen (EPR) experiment if it is possible to write the single and 
joint detection probabilities in the form |1J 

Pi = J p(\)Pi(\, <f>i)d\ ; p 12 = J p{\)Pi{\, </>i)P 2 (A, <h)d\, (1) 

A being the hidden variables with a probability distribution p(A). P»(A, 4>i) 
(i = 1,2), are some functions fulfilling the conditions < P»(A, <fii) < 1. <p\ 
and 02 represent controllable parameters of the experimental setup. 

Many experiments have been performed in order to test local realism vs. 
quantum mechanics via Bell's inequalities, but none of these experiments 
have violated Bell's genuine inequalities, which are based upon the assump- 
tions of realism and locality alone. All inequalities tested up to now involve 
additional assumptions, like no- enhancement. In the last few years [@ exper- 
iments using photon pairs produced in parametric down conversion (PDC) 
have become very popular. LHV models for EPR experiments with PDC 
have appeared in the literature ||, the aim of which has been to show that 
local realism had not been ruled out by these experiments, and to stress the 
relevance of the so called "loopholes" , in particular the one due to the low effi- 
ciency of optical photon counters. However, those models were mathematical 
constructs without physical content; they have no predictive power. 

The purpose of this paper is to exhibit a physical model which is based 
on the Wigner formulation of PDC. Our previous work Q has shown that 
the Wigner function of the electromagnetic field is positive for all performed 
experiments, and hence it provides a realistic description for the production 
and propagation of the radiation in these experiments. In particular, the 
Wigner function of the vacuum is the gaussian 

w(K,^)=n^- 2|ak|2 , (2) 

which in our local realistic model is interpreted as the probability distribution 
of a random radiation filling the whole space (the zeropoint field, ZPF). Thus 
the amplitudes {«k} play the role of the hidden variables A and W the role 
of the function p(A) in Eq.([L]). The nonlinear crystal, acted by a laser beam, 
gives an outgoing PDC field 
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E (+) = l-^s) [a^e-^ r+i ^ t + gale- l ^-^ r+i ^-^ t + -g 2 a^ r+ ^% 
k \ e o^o/ ' 2 

where the first term is the ZPF that crosses the crystal without any modifi- 
cation, and the other terms are produced via the non-linear coupling between 
the laser and the ZPF (g is the coupling parameter). 

The PDC beam propagates through a number of devices (lenses, beam 
splitters, etc.) characteristic of every experiment and finally some intensity 
lj arrives at every detector j. The intensity lj is a function bilinear of the 
amplitudes and o£ of Eq. (|Tj) . These functions lj are quite involved and 
have been studied in detail in Refs.0 for many experiments. The single and 
joint detection probabilities are given by 

Pj = J W({a k }, K})Q,(K}, K}, <f>)d N a k d N at , (3) 

Pia = J W({a k }, {«k})<5i({"k}, {a^}, 0i)<52({«k}, <p 2 )d N a k d N a£, 

(4) 

where 

Q = ^JdtJd 2 r [/({a k }, {a* k }, 0, r, t) - J ] , (5) 

1 = E + E~ being the intensity of the field arriving at the corresponding 
detector and 1$ is a constant which in our model is interpreted as the mean 
intensity of the the ZPF. The integration is carried over the time window 
and the surface aperture of the detector. We have divided [5] by the typical 
energy of one photon, so that Q becomes dimensionless, rj being the quantum 
efficiency of the detector. 

The relevant question is whether (§) and @ may be considered particular 
cases of (|1|) with Qj playing the role of Pj. The answer is not affirmative 
because we cannot guarantee that < Q < 1. Consequently we conclude 
that it is not possible to interpret directly the Wigner-f unction formalism as 
an LHV model for the PDC experiments. 

2 The Detection Model 

We shall devote the rest of the article to the description of the model of a 
detector that works in a strictly local way and provides the LHV model we 
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are searching for. We will first show the basic points of our model: 

1. The detector consists of a set of individual photodetector elements, D h 
each characterized by a frequency ui, and a wave vector k; (uji = |k;|/c), 
to which Di responds. We shall consider the direction of to be normal 
to the surface of the detector, which is taken as a cylinder of area ttR 2 
and length L. 

2. The relevant quantity for the detection is a filtered field corresponding 
to a detector element Df 

Hence it follows that the photodectector element D\ is sensitive to 
radiation with frequencies in the interval (07 — 4r, u)\ + 4p) with Alu « 
jt, T being the detection time window. If we assume that the incoming 
light beam has frequencies in the interval (^ m in,^max) with an average 
frequency uJ = {oJ max + u min )/2, and r is the coherence time of the 
beam, we shall have Sou = uj max — uj m i n ~ 2tt/t, so that the minimum 
number of detecting elements is N w 5uj/Auj w T/t. By putting typical 
values, T = 10~ 8 s and r = 10~ 12 s, we have the condition iV > 10 4 . 

3. We now define the effective intensity obtained from the filtered fields 
in the form 

JV 

7 = aoE^ +) ^f" ) . (7) 
1=1 

After that we replace the standard quantum Eq.(^) by the expression 

Q(7) = (l-e-« 7 - 7 °))e(7-/ m ) ; ( = V (hu)-\ (8) 

which completes the definition of our model. Iq is the average of / for 
the ZPF. I m is some threshold intensity fulfilling the condition I m > Iq, 
and Q(x) is the Heaviside function. 
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4. Now we rewrite the detection probabilities (^) and (f|) in the following 
equivalent form 

P = J P (T)Q(T)dI ; pi 2 = |pi 2 (7 1 ,7 2 )g 1 (7 1 )g 2 (7 2 )d7 1 d7 2 . (9) 

The probability distribution of the effective intensity is a gaussian, which 
is determined by the mean and the standard deviation (for details see ||). 
For instance, the probability distribution for the zeropoint field is 

The probability distribution of the effective intensity when there is a PDC 
signal present, p(I), is similar to that of the ZPF, but the following remarks 
are in order: (a) We shall define the signal mean effective intensity by I s = 
(I) — Iq. Here the word mean refers to both, time average over the detection 
window and ensemble average with the distribution function Eq.(0). If we 
assume that the signal enters parallel to the axis of the detector, it can be 
shown that the effective intensity, I s , is equal to the actual intensity, i.e. I s ~ 
I s . In practical situations the relation I s <C J is fulfilled. (6) The relation 
between the mean and the standard deviation is completely analogous to 
that of the ZPF alone (see Eq. fllOD ), and the distribution of the effective 
intensity I may be written 



\T T (I-I Q -Is) Z T 

v27r(i ) 

On the other hand, p(J 1 ,/ 2 ) is a doble gaussian function which is de- 
fined by the mean values of its marginals, their standard deviations and the 
correlation function ((7i - 7 ls - 7i )(7 2 - l 2s - 7 20 ))- 



3 The Detection Probabilities 

In this section we shall compare the predicted detection probabilities of our 
model with those of quantum optics. In the case of single detection proba- 
bility, let us consider the three following possible situations: 
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i) I s = 0. In sharp contrast with quantum optics, our model predicts the 
existence of some counts in any detector even in the absence of signal, the 
probability being very small if I m — I ^> a . 

ii) (I s <C o~o- This should be the normal situation in experimental prac- 
tice. In this case we may choose I m so that Io + I s — I m cr , but I m > I Q in 
order to preserve the positivity of Q in Eq. (|8|). With these two constraints 
it can be shown that p ~ (I s , a result that agrees with the quantum one if 
we put C _= i ■ 

in) (I s 3> o"o- In this case the detector saturates and gives a count in 
every time window. 

Finally, in the case of the joint detection probability the predictions of 
our model coincide with the quantum result in the limits ^-[1% — I%o) "C 1, 
which requires l is + 7 i0 - l im > &o- 

4 Constraints of the Model 

In our model there is a trade-off between the constraint Iq + I s — I m ^> a , 
required for the linearity of the response at low efficiency, and I m — 1 ^> cr , 
needed for the smallness of the dark counting probability. The need to satisfy 
both conditions implies that I s » o"o- This means that there is a minimal 
intensity of the signal which may be reliably detected, a constraint absent 
in the quantum theory of detection, but certainly existing in experimental 
practice. 

Let us analyze the consequences of the constraint. In experimental prac- 
tice a lens is placed in front of the detector in such a way that the signal field 
has spatial coherence on the surface of the lens. The condition for having 
spatial coherence is dX > RiRc, d being the typical distance between the 
nonlinear medium (with radius Rc) and the detector; Ri is the radius of 
the lens. On the other hand, the zeropoint field is not modified by the lens, 
which is evident because of the fact that energy cannot be extracted from the 
vacuum. As a consequence, the intensity of the incident signal is amplified 
by a factor b 2 = ir 2 Rf/\ 2 f 2 , f being the focal distance. On the other hand, 
84 % (91 %) of the total intensity is concentrated within the first (second) 
ring of the difraction pattern with a radius R = a x (f\/2Ri) = a\/A r , 
where A r = 2Ri/f is the relative aperture of the lens and a = 1.22 (2.23) for 
the first (second) ring . Consequently, the optimun radius of the detector is 
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given by R. 

By taking into account the above considerations, and using Eq. (10) and 
the condition dX > RiRc , the constraint I s ^> a gives 

R ate » Vf 2 R 2 c (12) 

which puts a lower bound on the single rates which may be used in reliable 
experiments. This result cannot be derived from (conventional) quantum 
theory. By putting typical parameters, that is 77 ~ 0.1, Ri, L and / of the 
order of fractions of a centimeter, 10 nm (which gives a coherence time r 
« 1 ps) we get a minimal counting rate of the order of 10 5 — 10 6 counts per 
second. This figure agrees fairly well with the actual experiments and should 
be considered a success of our model. 
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